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ABSTRACT 

In the collapsar model for common gamma-ray bursts, the formation of a centrifugally supported disk 
occurs during the first ~10 seconds following the collapse of the iron core in a massive star. This 
only occurs in a small fraction of massive stellar deaths, however, and requires unusual conditions. A 
much more frequent occurrence could be the death of a star that makes a black hole and a weak or 
absent outgoing shock, but in a progenitor that only has enough angular momentum in its outermost 
layers to make a disk. We consider several cases where this is likely to occur - blue supergiants with 
low mass loss rates, tidally-interacting binaries involving either helium stars or giant stars, and the 
collapse to a black hole of very massive pair-instability supernovae. These events have in common 
the accretion of a solar mass or so of material through a disk over a period much longer than the 
duration of a common gamma-ray burst. A broad range of powers is possible, 10 47 to 10 50 erg s , 
and this brightness could be enhanced by beaming. Such events were probably more frequent in the 
early universe where mass loss rates were lower. Indeed this could be one of the most common forms 
of gamma-ray transients in the universe and could be used to study first generation stars. Several 
events could be active in the sky at any one time. A recent example of this sort of event may have 
been the SWIFT transient Sw-1644+57. 

Subject headings: black holes - gamma-rays: bursts - stars:evolution, rotation - supernovae: general 



1. INTRODUCTION 

The collapsar model (Wooslcy 1993) was originally pro- 
posed as an explanation for common, "long-soft" gamma- 
ray bursts (GRBs) with durations near 20 s. This is a 
typical time scale for the gravitational collapse of the 
carbon-oxygen core of a massive star and the congru- 
ence of these time scales was an argument in favor of the 
model. It is well known, however, that forming a disk 
so quickly requires an unusually large amount of angular 
momentum to be s tored in the core up to t he point when 
the star dies (e.g.. IWooslev fe Hegerll2006ft . This is pre- 
sumably one reason why GRBs are rare compared with 
supernovae. 

Black holes are expected to fo rm in a significant frac- 
tion of massive star deaths (e.g.. lO'Connor fc Otdl201ll ) 
and it seems possible that, in at least some of these, 
the black hole will form without generating an outgoing 
shock. What observable signal would ac company such an 
event ? Would they all be "un- novae" (jKochanek et al.l 
2008), stars that simply suddenly disappear from view? 
Or will some sort of transient display announce the birth 
of any black hole? Here we explore the possibility that 
these black hole systems form an accretion disk, but only 
after almost all of the star has accreted into the hole. 
Torques in differentially rotating contracting stars tend 
to concentrate angular momentum in the star's outer lay- 
ers. The surface will also have a large angular momentum 
in a tidally locked system. On the other hand, stellar 
mass loss preferentially removes these same layers, so, 
just as in the usual collapsar model, reducing mass loss, 
e.g., by reducing the metallicity, favors the retention of 
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high angular momentum at the stellar surface. 

Because the outer layers of the star take longer to fall 
in - months to years in the case of a red supergiant - the 
duration of the transient produced by accretion through 
a disk in these cases is much longer than for the common 
GRBs studied so far. They might even be confused with 
other forms of gamma-ray transients like those arising 
from super-mas s ive black hole accreti o n in active galaxies 
(IWooslevI 120111: iBurrows et alj I20H iLevan et all 120111: 
iQuataert fc Kasenll2011l b 

Here we shall consider four varieties of long duration 
transients categorized by the kind of system where they 
occur: 1) single supergiants with metallicity less than 
about 0.1 Zq and low mass loss (§ 13. 1|) : 2) red or blue 
supergiants in tidally-locked, or nearly tidally-locked bi- 
naries (§ I3.2[) ; 3) pair- instability induced collapse to a 
black hole (§ 13 . 3|) : and 4) helium stars in tidally locked 
binaries (§ 13. 4[) . In each case, the duration of the tran- 
sient that is produced is correlated with the radius of the 
presupernova star, hence shortest for the helium stars, 
while the luminosity is inversely correlated with the ra- 
dius, hence faintest for the red supergiants. As with 
common GRBs, we are unable to predict ab initio the 
spectrum of these transients or the efficiency for convert- 
ing black hole accretion rate into electromagnetic radia- 
tion. Presumably, as in other forms of transients powered 
by very rapidly accreting, rotating black holes, accretion 
energy and black hole spin energy power relativistic bi- 
polar outflows. Poynting flux, shocks within these out- 
flows, and shocks with the circumstellar medium make 
x-rays and gamma-rays. It is common to take an effi- 
ciency of 1 % to 10% times Mc 2 for the efficiency (e.g., 
McKinney 2006) thus suggesting an integrated power of 
~ 10 52 erg — 10 53 erg per solar mass accreted. Interest- 
ingly these "failed supernovae" (jWooslevI I1993D are far 
more brilliant and powerful than their (model-wise) more 
successful cousins, ordinary supernovae. 
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The idea that late time accretion in supernovae that 
make black holes might power long ga mma-ray (or hard- 
x-ray) transients was suggested by iMacFadven et al.l 
(|2001h . The transients discussed there, referred to as 
"Type 2 collapsars" , were produced because of fall-back 
in an explosion that launched a relatively strong outgoing 
shock in a red supergiant progenitor. Two dimensional 
models explored jet propagation and break out for a vari- 
ety of fall-back rates, but little attention was given to the 
angular momentum distribution necessary to form a disk. 
Here we consider a wider variety of progenitors and pay 
closer attention to t he ang ular momentum, especially in 
the outer layers. Q ()2003l ) also discussed the possibility 
of very-long duration X-ray transients powered by fall 
back and accretion through a disk. He suggested these 
might look like "ultra-luminous x-ray sources" seen in 
nearby galaxies. These would be much longer in dura- 
tion (thousands of years) and have much lower luminosity 
than the events considered here, but are an example of 
the same genera l idea. 

Mor e recently iQuataert fc Kaseiil (|2011l) and IWooslevI 
(|2011h suggested that fallback or a failed explosion in 
a red supergiant could give a long duration gamma- 
ray transient that might be confused with a gamma-ray 
blazar. Quataert and Kasen, in particular, suggest this 
as an alt ernate interpretation of SWIFT transient Sw- 
1644+57 (|Burrows et alllMtl . Woosley's paper was ac- 
cepted prior to the discovery of this event. 

A variety of other models may lead to long dura- 
tion bursts. For example a white dwarf (jFrver et al.l 
1999) or helium star merging with a black hole or neu - 
tron star dFrver fe Wooslevl Il998t IZhang fc Frverl 120011) . 
A neutron star or black hole merging with a a giant 
star might give a Thorne -Zytkow object (|Terman et al.l 
[i~995l lPodsiadlowskill2007l) or a black hole accreting at a 
rapid rate (jChevalierl 1 1 993f ) . Helium star merger with a 
compact object also gives a long transient of some sort 
(jFrver fc H cgcr 2005|). Here we consider only those tran- 
sients coming from stars that are not actively merging 
with a companion. 

2. THE PRODUCTION OF STELLAR MASS BLACK HOLES 
AND TYPE 3 COLLAPSARS 

A black hole will form inside a massive star when a 
sufficient concentration of matter and energy accumu- 
lates to form an event horizon. The simplest possibil- 
ity is that some portion of a massive star experiences 
an instability and falls freely, or almost freely, until it 
reaches such a high density that a trapped surface is 
formed. This rarely happens in nature though. Only 
in t he most massive stars e xperiencing a pair instabil- 
ity (IHeger fc Wooslevl 2002 ) or a general re lativistic in- 
stabilitv (jFuller et al.l fl986t iMontero et all 1201 ll ) can a 



trapped surface form before the center reaches nuclear 
density and bounces. That is, black hole formation in a 
massive star is almost always preceded by protoneutron 
star formation (jO'Connor fc Otdl20ll . 

For stars in the range ~ 20 M Q — 100 Mq, however, 
the protoneutron star that is formed experiences rapid 
accretion from the surrounding layers of heavy elements 
- silicon, oxygen, etc. If rotation, magnetic fields and 
neutrinos do not readily produce an outgoing shock, the 
hot protoneutron star grows beyond a critical mass and 
collapses to a black hole. The accretion rate is larger 



for more massive helium cores (|Frverl 11999ft . and black 
hole formation without an outgoing shock may be a fre- 
quent occurrence in stars with main sequence masse s 
over about 20 M (jO^Co nnor fc Qttl[201lt lSmartdl2009h . 
Whether all stars above about 20 M produce black 
holes without outgoing shock waves is uncertain, but the 
possibility is not precluded by observations or theory. 
"Prompt" black hole formation is probably more likely 
in very low metallicity stars that have experienced little 
mass loss and thus have especially large helium cores, 
but it could occur for solar metalli city stars, depend ing 
upon the uncertain mass loss rate (|Smith et al.ll201ll ). 

It is also possible that such a weak outgoing shock 
is generated that most of the star fails to explode an d 
a black hole is made by fall back (IZhang et al.l [2 008) . 
These were the systems explored by IMacFadven et al.l 
(|200l . The ones treated here differ in assuming that any 
outgoing shock is insufficient to eject more than a small 
fraction of the surface layers of the star. We shall call 
such events, where only the surface layers have enough 
rotation to form a disk and a shock so weak that these 
layers are not ejected "Type 3 Collapsars". This is to 
distinguish them from ordinary GRBs where the disk 
forms promptly from material deep insi de a highly dif- 
ferentially rotating helium co re (Type 1: IWooslevI Il99l 
IMacFadven fc Wooslevl I1999D and from systems where 
longer duration transients result from the delayed fall 
back of other material also deep inside th e star (Type 2; 
IMacFadven et~aT]|200ll: IHeger et al.ll200l . 

3. EXAMPLES OF TYPE 3 COLLAPSARS 

The preservation of angular momentum in the surface 
layers of the star can occur either because mass loss is 
inefficient in a single star or because of tidal interaction 
with a close companion. Here we consider examples of 
both sorts. 

3.1. Massive Stars with Low Metallicity and Low Mass 

Loss 

A paucity of metals in a star affects its evolution in sev- 
eral ways - reduced energy generation by the CNO cycle 
in the hydrogen burning shell, reduced opacity, and re- 
duced mass loss. Of these effects, the reduction in mass 
loss is least well understood, but the lack of atomic lines 
in the atmosphere on the main sequence and of grains 
as a red giant seems certain to lead to some reduction. 
A common assumption is that the mass loss for mas- 
sive main sequence stars depends upon s ome fractional 
power of the metallicity (Kudritzki 2002). More impor- 
tant and less certain is the dependence of mass lost as a 
red or blue supergiant on metallicity. It is thought that 
mass loss from cool giants is more dependent upon pul- 
satio ns and grain formation (|Reimcrs 1977; Smit h~et al.l 
1201 lh . One might expect therefore a rapid fall off in mass 
loss below some value necessary for significant grain pro- 
duction. It has been estimated that r ed giant mass loss 
will be significantly less below 0.1 Z Q (|Bowen fc Willsonl 

ll22llz0kk3l2QQ4[)- 

If stars do not lose mass then they conserve angular 
momentum. Since the natural course of evolution leads 
to the contraction and spin up of the inner star, shear 
instabilities and magnetic torques will concentrate an 
increasing amount of angular momentum in the outer 
regions of the star. Whereas the inner part does spin 
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faster due to contraction, it actually loses most of its 
initial angular m omentum by the time it r eaches central 
carbon burning. iHeger fc Wooslevl (|2010f ) surveyed the 
evolution of non-rotating zero-metal massive stars from 
10 M^ to 100 M(7). In a p arallel study currently underway 
(|Heger fc Wooslevi[201lll . we are surveying the evolution 
of rotating stars in the same mass range with metallic- 
ity 0, 10~ 3 and 10 _1 that of the sun. The two stars 
discussed here are extracted from that survey and are 
typical in the way that they accumulate large amounts 
of angular momentum in their outer layers before dying 
as supergiants. All calculations of s tellar evolution pre- 
sented here used the Kepler code (IWeaver et al.l 119781 : 
IWooslev etai1l2002l;lHeger et al.ll2000D and includedan- 
gular momentum transport by magnetic torques (jSpruitl 
l2002t ). Models V24 and V36 are 24 M Q and 36 M main 
sequence stars with metallicity 0.1% solar. Both stars 
rotated rigidly on the main sequence with a moder- 
ate speed of about 200 kms -1 which is about 20% of 
the Keplerian value. Both presupernova stars had hy- 
drogenic envelopes that were, throughout most of their 
mass, radiative. Model V36, however, had a low mass 
surface convection zone that included 0.11 M . Model 
V36 was consequently a yellow supergiant at death (L — 
2.2 x lO^ergs" 1 ; R = 5.4 x 10 13 cm; T cS = 5700 K). 
Model V24 was a rather large blue supergiant (L = 
1.1 x 10 39 ergs-\ R = 1.0 x 10 13 cm; T cS = 11, 300 K). 

These two stars ended their lives with cores of he- 
lium and heavy elements of 8.3 M and 15.0 M Q respec- 
tively and angular momentum distributions as shown in 
Fig.[TJ Whereas the lack of sufficient angular momentum 
within the helium core precludes making a disk around 
a black hole within that mass, there is ample rotation in 
the outer part of the hydrogen envelope to do so. This 
could be a very common occurrence. From the survey of 
IHeger fe Wooslevl t|201 lh . the total angular momentum 
of a massive star at birth that has an equatorial rotation 
speed of about 20 % Keplerian on the main sequence is 
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If the star does not lose mass, most of this angular mo- 
mentum becomes concentrated, in the presupernova star, 
in a nearly rigidly rotating hydrogenic envelope with ra- 
dius either < 10 13 cm (blue supergiant) or ~ 10 14 cm (red 
supergiant). Whereas the density declines with radius 
in the actual envelope, one can obtain some interesting 
scaling relations by assuming constant density. For a mo- 
ment of inertia, / w 0.4Mi? 2 , the angular velocity at the 
surface, R — 10 14 Ru cm, is 



1.2 x 10" 10 i?; 



M 



20 M m 



0.8 



rads 



(2) 



The specific angular momentum of this rigidly rotating 
envelope is j = ^ojr 2 , or about 10 18 cm 2 s _1 at the sur- 
face of a red supergiant and 10 20 cm 2 s -1 for a blue su- 
pergiant, if no mass is lost. 

These can be compared with the angular momentum 
required to make a disk around a non-rotating black hole, 
j = 2V3GM/c = 3.1 x lO 17 (M BH /2OM )cm 2 s- 1 or 
j = 2/VZGM/c = 1.0 x 10 17 (M B H/20M Q )cm 2 s- 1 for 
a maximally rotating hole (Kerr parameter a = 1). All 



massive stars that do not lose mass have sufficient angu- 
lar momentum in their outermost layers to make a disk 
around any black hole formed in their collapse. 

Of course stars do lose mass and the fact that, for 
constant u), j oc r 2 means that the first mass to be lost 
contains most of the angular momentum. In practice, 
we find that a red supergiant that loses more than a few 
percent of its total mass before dying will probably not 
make a black hole-disk system. 

This sensitivity to mass loss is illustrated by the evo- 
lution of a 25 Mq model, 025, with 10 % solar metallic- 
ity in a calculation that included mass loss at one-half 
the standard value for that metallicity. The star lost 
1.17 M (5%) of its mass and died as a red supergiant 
with a radius of 1.1 x 10 14 cm. The value of j in its out- 
ermost layers was 2.8 x 10 17 cm 2 s^ 1 , too little to form a 
disk around a Schwarzschild black hole, but possibly just 
enough to form a disk around a mildly rotating one. The 
point where j declined to one-half its surface value in 
the presupernova star was located 5.5 M interior. This 
is borderline case. Stars with less mass loss, i.e., lower 
metallicity, and stars that never become red giants could 
make long gamma-ray transients (Table [TJ . 



3.2. Supergiants in Interacting Binaries 

Based upon theoretical models for single massive stars, 
the red supergiant progenitors of common Type Hp su- 
pernovae are not expected to have rapidly rotating en- 
velopes. Expansion to a red giant slows the surface lay- 
ers. Torques transfer the angular momentum of the more 
rapidly rotating inner regions outwards concentrating it 
in the outermost layers. For stars of sufficient metallic- 
ity, stellar winds are efficient at removing these layers 
and, along with them, most of the star's angular mo- 
mentum. For example, a 25 M solar metallicity star 
rotating rigidly with an equatorial speed of 200 kms -1 
on the main sequence ends up a 12.1 M supernova pro- 
genitor if commonly used estimat es for mass loss are 
employed (jWooslev fc Hegerl l2007h . The surface rota- 
tion rate of that red supergiant (R = 9 x 10 13 cm) is 
2 x 10~ 12 rads -1 and its angular momentum at the sur- 
face is 1.6 x 10 16 cm 2 s _1 . Nowhere in the star does 
any matter have enough angular momentum to form a 
disk around a black hole. If a black hole forms without 
making an outgoing shock, the star abruptly disappears. 
This sort of "un-nova" might be a common o ccurrence 
in the modern universe (jKochanek et all [2008). though 
see lSmith et~aT1 (|2011[ ). 

The outcome can be quite different in a binary sys- 
tem, however. A s tellar merger could lead to a ra pidly 
rotating envelope (jlvanova fc Podsiadlowskil l2003f). but 
such mergers probably lead to ejection of a common en- 
velope. An interesting alternative is detached binary sys- 
tems with a small enough separation for tidal interaction, 
but too distant to merge or exchange mass. 

The circularization and synchronization time scales 
for binary stars, one of which is a red giant, have 
been considered by IZahnl (|1966l ). IZahnl (|1989l ) and 
IVerbunt fc Phinnevl (|1995f) . For a primary with mass, 
M, envelope mass M env , luminosity L, and radius R, 
in an orbit having radius d, with a companion of mass 
M2, the circularization time, r c i rc is (|Verbunt fc Phinnevl 
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where / is a number close to unity. The lifetime of a 
25 M© main sequence star as a red supergiant is approxi- 
mately 700, 000 yr, depending upon uncertain convection 
parameters. The mass of the envelope that is convective 
at any one time, however, is only about 3 M Q and the 
star loses more than this in its lifetime. A more appro- 
priate time scale might be the time the star spends losing 
its last 3 Mq, which is about 250, 000 yr. The mass of the 
star when it becomes a supernova is 12.1 M© and we as- 
sume a companion mass of, e.g., 6M©. The luminosity 
of the star as a red supergiant is 1.0 x 10 39 ergs _1 and 
its moment of inertia is approximately 4 x 10 61 gcm~ 2 . 
Equation ([3|) then implies that during its last 250, 000 yr, 
orbits with radius d < 3.8R s» 18 AU will circularize. 

The synchro nization timescal e for co-rotation is some- 
what different (|Zahnlll 9771 . 1 19891) and, in the present case, 
is easier to achieve. 
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Most of the angular momentum is in the hydrogen en- 
velope and, from the same 25 Mq model, MR 2 /I = 2.9. 
With q ss 1 and A sy 



0.02 (lZahnlll989h . 
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which only requires d < 6.7i? for synchronization in 
250, 000 yr. In such a system, the period will be shorter 
than about 30 years (w = 10 _8 rads _1 ). Systems much 
closer will merge or experience mass exchange, so we take 
this to be a typical rotation period for the envelope in 
a detached binary system that might achieve corotation. 
Longer period systems will still transfer significant rota- 
tion to the envelope, however. 

This is about four orders of magnitude more angular 
momentum in the outer envelope than in the case of the 
single star. Since the outer several solar masses of the 
presupernova star is convective, this rotation rate will 
persist to some depth. All of that mass would have suf- 
ficient angular momentum to form an accretion disk if 
its helium core collapsed to a black hole. Even then 
though, the total angular momentum of the envelope, 
2.0 x 10 53 ergs, would be a small fraction of the orbital 
angular momentum of the binary pair of stars. The addi- 
tion of only 1 % of the angular momentum estimated here 
would still result in large quantities of material forming 
a disk, so the range of binary separations allowed is ac- 
tually larger, perhaps by a factor ~ (100) 1 / 6 , or about 
two. 

Since the time scale for material to fall in from 10 13 — 
10 cm is longer than any realistic viscous time scale for 
the disk, the accretion rate will be given by the collapse 
timescale of the envelope. For matter at the base of the 
hydrogen shell where a disk would first form, the en- 
closed mass is 9.3 M© and the radius about 1 x 10 13 cm. 
At the outer edge of the star the mass is 12.1 M© and 
the radius, 9.8 x 10 13 cm. The corresponding free fall 

times (r ff = (2A<KGp)~ 1/2 = 446 (p/1 gcuT 3 )~ 1/2 s) are 



2.1 x 10 5 s and 5.7 x 10 6 s suggesting that the event would 
last for months. Accreting 3 M© in 10 6 s with 10 % effi- 
ciency for converting accreted mass to outgoing energy 
would give a jet power of 10 48 erg s -1 . Conversion of only 
10% of the jet power (1 % of the accreted mass energy) 
into gamma-rays would explain the recently discovered 
transient Sw 1644+57 (Bu rrows et al.ll2011l) . especially 
if moderate beaming were involved. 

Whereas the focus here is on red supergiants because 
they give the longest transients, tidal locking would also 
occur in blue supergiants with smaller radii. The key 
point is that forced co-rotation with a body in Keplerian 
orbit gives, for a large range of orbital separations, enve- 
lope material that is rotating too fast to accrete, with- 
out hindrance, into a stellar mass black hole. Since the 
time scale for accretion is set by the collapse time for the 
outer part of the star, not the viscous time of the disk, a 
transient from a tidally locked blue supergiant would be 
similar in duration within a factor of a few to those stud- 
ied for single blue supergiants in § 13.11 Blue supergiants, 
however, will have less mass in the outer layers of the 
envelope and hen ce should have less powe rful outbursts 
(Appendix [A] and iQuataert k, Kasenll2011l ). 

3.3. Pair-Instability in Very High Mass Stars 

Following helium burning, non-rotating helium cores 
with masses greater than 133 M Q will collapse di- 
rectly into black hol es with no outgoing shock 
(|Heger fc Wooslevll2002[ ). The collapse is caused by the 
pair-instability and, above this mass, nuclear burning is 
unable to reverse the implosion. Rotation, unless it is 
very rapid and highly differential, does not change this 
mass limit greatly since the instability exists in the deep 
interior where the ratio of centrifugal force to gravity is 
small. 

The collapse of a rotating pair unstabl e star to a 
black hole was followed in two dimensions by lFrver et al.l 
(2001). First, a 300 M© main sequence star was evolved 
to the point of instability in the KEPLER (ID) code. 
This star produced a helium core of 180 Mq. The initial 
star was assumed to rotate rigidly on the main sequence 
with a ratio of equatorial speed to Keplerian of 20%. An- 
gular momentum transport, especially by shear mixing 
and Eddington Sweet circulation, was followed through- 
out the evolution, but magnetic torques were neglected. 
As a result, the collapsing helium core rotated sufficiently 
rapidly to produce a disk outside of a massive black hole 
once about 140 M© of the 18 Mq core of heliu m and 
heavy elements had collapsed. iFrver et all (|2001h specu- 
lated that the accretion of the remaining 30 M© — 40 M Q 
might produce a very bright, lengthy GRB. 

We have recalculated th e evolution of similar stars, 
but with magnetic torques (|Spruitll2002l ) included in the 
simulation. Models Z250A, B, and C are stars with 
zero metallicity and a mass at the time they formed 
of 250 M©. They differ in the amount of angular mo- 
mentum each had on the main sequence. Z250A had 
0.75 x 10 54 ergs corresponding to a rotation rate of 
170 km s _1 , or 12 % Keplerian. Models Z250B and Z250C 
had 1.0 x 10 54 erg s and 1.5 x 10 54 erg s of angular momen- 
tum and rotational speeds of 220 kins" 1 (15 % Keplerian) 
and 310 km s -1 (19% Keplerian). These made helium 
cores of 143 M Q , 166 M Q , and 222 M© respectively, all of 
which collapsed directly to black holes following helium 
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depletion. Unlike in lFrver et al.l ()2001l ) however, none of 
these helium cores had sufficient angular momentum to 
make a disk around a black hole anywhere in their inte- 
rior. Because mass loss was neglected in these models, 
the angular momentum lost from the helium core due to 
magnetic coupling instead ended up concentrated in the 
outer part of the hydrogen envelope. 

Fig. [3] shows the evolution of the angular momentum 
in Model Z250B. The evolution of the other two models 
was similar. Inadequate angular momentum exists in the 
helium core to form a disk but, for Model Z250B, the 
outer 16 M can still do so. As in the tidally locked red 
supergiant case, the collapse time scale for this material 
is days to months, so a high energy transient of around 
10 47 ergs -1 , possibly boosted by beaming, could exist for 
that time. 

It should be noted, however, that mass loss is even 
more uncertain in these models than in the others. Due 
to rotationally induced mixing, the envelopes of all three 
models were appreciably enhanced in C, N, and O and 
thus were red supergiants when they died. Neglecting 
mass loss in such a situation is probably not justified. 
Also there may be other mechanisms for losing mass in 
such massive stars besid es line and grain driven mass loss 
(|Smith k Owockill2Ql06ft . 

3.4. Helium stars in Binaries 

Helium stars in close tidally locked binaries have been 
frequently suggested as progenito r s for common GRBs 
(iTutukov k Cherepashc huk 2003; Podsiadlows ki et all 
20041: llzzard et all 12004 iBogomazov et al.l 120071: 
van den Heuvel k Yoonl |2007|) . On the positive side, 
such events should occur frequently in nature and, if 
the core of a massive Wolf-Rayet star collapses to a 
black hole while the surface is forced to rotate at a 
fraction of its Keplerian speed, disk formation around 
a black hole is assured. On the negative side, however, 
the production of a GRB with typical duration ~10 s, 
requires that the inner part of the helium star rotate 
rapidly, not just its surface. The inclusion of magnetic 
torques in tidally locked stars leads, by itself, to too 
little rotation to make a disk in that part of the star 
which might accrete rapidly ()Wooslev k Hegerl 120061 : 
Ivan den Heuvel k Yoonll2007t ) . On the other hand, leav- 
ing out magnetic torques makes it difficult to produc e 
the spin rates of ordinary pulsars (He ger et al.l [2005). 
and could lead to an overabundance of GRB progenitors. 

This is not to say that close binary systems do not pro- 
duce some form of gamma-ray transient though. The an- 
gular momentum in their cores may be adequate to pro- 
duce a millisecond magnetar, and thus a common GRB, 
just not a (Type 1) collapsar. If a black hole forms, the 
accretion of the outer layers of the star may lead to the 
formation of a disk and GRB much longer than typical. 
It is this lat ter possibility that we explore here. 

Following Ivan den Heuvel k Yoonl 12001 . we consider 
the evolution of helium cores of 8 M and 16 M whose 
surface s have bec ome tidally locked with a clos ely or- 
biting companion. Ivan den Heuvel k Yoonl IpOOl found 
the shortest possible periods were 2.05 hours (uj — 8.5 x 
10- 4 rads~ 1 ) and 2.47 hours (lu = 7.1 x lO^rads- 1 ) 
for 8M and 16 M respectively, and they gave exam- 
ples of systems containing a compact object and a helium 
core that might approach these limiting values. An 8 M 



He core in a tidally locked binary with a 0.8 M com- 
panion filling its Roche lobe had a longer period of 7.17 
hours (u> — 2.4 x 10~ 4 rads _1 ). Here we calculate two 
models for each helium core mass. Model 8A has an en- 
forced rotation rate in its outer 0.5 M of 2 x 10 -4 rad s _1 
throughout its evolution. Models 8B, 16A, and 16B have 
rotation rates of 8.5 x 10 _4 rads _1 , 2 x 10~ 4 rads _1 , and 
7 x 10 _4 rads _1 respectively. 

Fig. |4] shows the evolution of the angular momentum 
distribution as a function of mass for these four models. 
During helium burning, magnetic coupling maintains a 
state of nearly rigid rotation with an angular speed given 
by the surface boundary condition. For both Model 8A 
and 16A, this corresponds to a ratio of rotational speed 
to Keplerian of ~ 8%. For Models 8B and 16B, the ra- 
tio is - 32%. Roughly the inner 75% (8M ) or 85% 
(16 Mq) of the helium star is convective during helium 
burning and rigid rotation is enforced in that region. The 
outer part of the star is radiative, but magnetic torques 
there suppress differential rotation. During carbon burn- 
ing and more advanced stages, the inner core contracts 
and rotates more rapidly. Magnetic torques are too inef- 
ficient to enforce rigid rotation. Later, during carbon 
burning, the stars ignite a convective helium burning 
shell which persists until the star dies. This convective 
shell does not extend all the way to the surface though, 
but does extend beyond the 0.5 M at the surface where 
corotation is enforced. Thus all stars die with an an- 
gular velocity in their helium convective shell given by 
the surface boundary condition, but have larger rotation 
rates in their cores, due to contraction during advanced 
burning stages. 

Fig. [5] shows the corresponding angular momentum in 
the evolving models. Shear instabilities and magnetic 
torques transport angular momentum out of the inner 
core and to the surface. In the absence of a binary com- 
panion, the surface would have actually spun up, but 
would have been removed by mass loss. In the present 
parametrized calculations though, there is no mass loss 
and the set rotation rate of the surface absorbs any excess 
angular momentum delivered to it. Consequently the to- 
tal angular momentum of the star decreases with time. 
For Model 8A at helium ignition, it was 9.0 x 10 51 ergs, 
but for the presupernova star it was 1.6 x 10 50 ergs. 

Fig. [5] shows the final distribution of angular momen- 
tum in the stars compared with that required to make a 
disk around a black hole, once the mass interior to the 
given sample point has collapsed. Four curves are given. 
Three are for the minimum angular momentum required 
to form a disk at the last stable orbit of a) a non-rotating 
(Schwarzschild) black hole; b) a Kerr black hole; and c) 
a black hole with rotation given by the actual angular 
momentum distribution of the presupernova star. The 
fourth curve is the angular momentum distribution in 
the presupernova star. Because the inner parts of the 
star are not rotating very rapidly, curve a) and c) are 
similar except near the surface. Only where the actual 
angular momentum lies above curve c) can a disk form. 
Material with less angular momentum plunges directly 
into the hole. 

None of the models has enough angular momentum to 
form a disk in its inner core. Thus a Type 1 collapsar will 
not occur. The angular momentum in the inner 1.7 M 
of the stars is still appreciable, however. Models 8A and 
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16 A would form 1.4 M Q (gravitational mass) pulsars with 
periods of 7.4 ms and 2.6 ms respectively. Models 8B and 
16B would form pulsars of 3.2 ms and 1.5 ms. Model 16B 
thus qualifies as a "millisecond magnetar" candidate and 
might power a common GRB if a pulsar is able to form 
before accretion turns it into a black hole. Models 8B, 
16A and 16B could also make powerful pulsar-powered 
supernovae. Model 8A and other more slowly rotating 
systems probably would not. 

Of greater interest for the present paper, however, 
is the large angular momentum in the outer layers of 
three of the models. These would be relevant if the 
center of the star collapsed to a black hole rather than 
a millisecond pulsar. Model 8A has angular momen- 
tum sufficient to form a disk in its outer 0.074 M Q . 
Two other models (Table [T]) have similar large rota- 
tion in their surface layers. This is cons istent with 
the results of Ivan den Heuvel fc Yoonl (|2007l) who found 
that their most rapidly rotating 8 M Q and 16 M Q he- 
lium cores ended their lives with j in outer layers of 
3.7 x 10 17 cm 2 s~ 1 and 6.0 x 10 17 cm 2 s _1 . The radius 
and gravitational potential of the high angular momen- 
tum surface layers allows an order of magnitude estimate 
of the free fall time scale (Table [1]), about 100 s. 

The tidal interaction actually acts to brake the rota- 
tion of the surface layers compared to what they would 
have had in a star without mass loss and no compan- 
ion star. To illustrate this, Model 8C was calculated, 
starting from the rigidly rotating helium burning stage 
of Model 8B (total angular momentum 3.6 x 10 51 ergs), 
but with no surface boundary condition and no mass 
loss. That is the rotation rate of the surface layers was 
allowed to adjust to be consistent with whatever angu- 
lar momentum was transported to them. Without mass 
loss, the star was forced to conserve angular momentum 
overall and the very outer layers ended up rotating very 
rapidly (Fig. [6]). The angular speed, rather than being 
u = 7.5 x 10~ 4 rads _1 was 3.7 x lO^rads^ 1 . In fact 
the outer tenth of a solar mass rotated so rapidly that 
it would be centrifugally ejected. In a realistic calcu- 
lation, this matter and more underlying matter would 
probably have been ejected as a disk or a centrifugally 
boosted wind. The deep interior still lacked sufficient an- 
gular momentum to form a disk around a black hole, but 
the amount of surface material that could form a disk 
(if mass loss did not remove it) was significant, about 2 
Mq, comparable to that in Model 8B. Two other models 
like 8C were calculated that included mass loss appro- 
priate for a Wolf-Rayet star at solar metallicity and one- 
tenth that value. The first ended up with a final mass 
of 3.95 M Q and no surface layers with sufficient angular 
momentum to form a disk. The second ended up with a 
mass of 7.09 M Q and still had enough angular momentum 
in its outer solar mass to make a disk. 

So the major role of the close companion in the calcu- 
lations done here is to preserve the angular momentum 
in the outer layers that might otherwise be removed by 
mass loss. A companion may also be essential so that 
the helium core rotates rapidly in the first place, fol- 
lowing the loss of the star's hydrogen envelope. When 
magnetic torques are included in the calculation, the he- 
lium cores of supernova progenitors that are red and blue 
supergiants rotate too slowly for any portion of that core 
to make a disk. A close companion could either preserve 



angular momentum in the core by removing the envelope 
very early in the evolution or add angular momentum by 
tidal locking later on. 

4. CONCLUSIONS 

Table [1] and Table [2] summarize the wide variety of 
transients that can be produced by Type 3 collapsars. In 
Table [1] Afdisk is the amount of mass capable of forming 
a centrifugally supported disk around a black hole with 
a Kerr parameter given by accreting all the mass interior 
to Mdisk and R is the range of stellar radii where that 
mass exists. The accretion time scale has been estimated 
using the radius for the center of mass of -M^sk- The lu- 
minosity has been estimated by multiplying the effective 
accretion rate (mass/time scale) by 1 % c 2 . This assumes 
10 % efficiency for making a jet and 10 % efficiency for 
converting jet energy into hard radiation. Most mod- 
els produce extreme Kerr black holes (a = 1, Tabled, 
so th e radiative efficiencies may be higher than assumed 
here (|Tchekhovskoy et al.l 12011). The numbers in the 
table are only order of magnitude estimates. 

Table [2] gives the the estimated pulsar period formed 
in the stellar collapse assuming the angular momentum 
that exists in the inner 1.7M of the core, a final grav- 
itational mas s of 1.4 M(T), and a moment of inertia of 
0.35M grav i? 2 (|Heeer et al.ll2005D . If a black hole forms 
before the protoneutron star has radiated away most of 
its gravitational binding energy, the rotation rate will 
be slower as will the rotational kinetic energy available 
for explosion. Af prc -sN and Mno-core are the masses of 
the presupernova star and its helium core respectively 
and M(j > j cr it) is the Lagrangian mass coordinates of 
the stellar layers that can form a centrifugally supported 
disk. The Kerr parameter, a, is what would develop if the 
entire presupernova, including its rapidly rotating outer 
layers, collapsed to a black hole. In case where it exceeds 
unity it has been capped at 1. The Kerr parameter for 
the helium core alone is an c . BE(j > j C rit) is the net 
binding energy, in units of 10 49 erg, of the material that 
could form a disk. 

A range of masses can form disks and a large range of 
time scales characterize the accretion, hence the possi- 
ble luminosities span many orders of magnitude. The 
shorter brighter transients from Wolf-Rayet stars in 
tidally locked binaries might comprise an interesting sub- 
set of long, but otherwise ordinary GRBs, as others have 
noted. This channel could be especially important in 
regions where high metallicity and high mass loss pre- 
clude the production of GRBs by solitary stars. The 
mass of 56 Ni produced in the collapsar model is sensi- 
tive to the temperature, density, time scale, and ejec- 
tion fraction of the disk and is difficult to estimate. The 
lower accretion rate, larger black hole mass, and smaller 
Kerr parameter for the models presented here, however, 
will give lower temp eratures and densities in the disk 
(|Popham et al.lll999fh Below about 0.01 M s" 1 , the na- 
ture of the accretion changes and the disk is no longer 
neutrino dominated. A thicker disk might qualitatively 
change the nature of the transient (MacFadyen, private 
communication). It is thus reasonable to expect that 
the 56 Ni production, per solar mass of material accreted, 
will be different here and perhaps much smaller. These 
bursts might not be accompanied by supernovae (e.g., 
iGal-Yam et aLll200l . 
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On the longer end, blue and red supergiants will pro- 
duce transients of order 10 4 s and 10 5 s respectively. 
These could result from either tidally locked systems 
or more distant systems with low metallicity and mass 
loss. These fainter and potentially more common systems 
may not have been detected yet and certainly should be 
searched for. They could easily be confused with flaring 
gamma-ray blazars (lWooslevll2011l; iBurrows et al.l l2011t 
iQuataert fc Kasenl 120111 : iLevan et al. I12011D, but would 
differ in that they only happen once and would be asso- 
ciated with regions of massive star formation, not par- 
ticularly galactic nuclei. They too would probably not 
be accompanied by supernovae and might have a differ- 
ent beaming factor and spectrum from ordinary GRBs. 
The black holes produced in these supergiant progenitors 
would frequently be rotating at near the maximum value 
(Fig. [1] and Fig. HD. 

An exciting prospect is that the first generations of 
stars, those with metallicities from zero to 0.1 solar might 
be studied directly using the hard x-rays they emit when 
they die. Whereas predictions are still clouded by the 
vagaries of mass loss and the explosion mechanism, our 
studies predict that the fraction of massive stars emitting 
such bursts could be substantial, essentially limited only 
by the fraction that make black holes promptly. If none 
make black holes promptly, then the death of every mas- 
sive star must make a supernova, though perhaps a faint 
one. This too is an interesting possibility worth further 
exploration. 

Because of their long duration and potentially high 
event rate, it is possible that several Type 3 collap- 
sars will be yisible in distant galaxies at any one time. 
iMadau e~ al. (1998) estimate 20 core-collapse supernovae 
occur each year in a field of view 4' x 4' in the redshift 
range z = Q . . . 4. This is about 6 supernovae per second 
spread over the whole sky. The observation of GRBs 
beyond redshift 8 shows that massive stars were dying 
at an even earlier time. Black hole formation proba- 
bly occurs in a significant fraction o f these supernovae 
(jO'Connor fc Ottlboill: lSmarttl [2009f). perhaps all those 
with mass over 20 M Q and metallicity less than 10 % so- 
lar. In order to estimate, crudely, an event rate, as- 
sume that prompt black hole formation occurs in 10 % 
of core collapses and 1 % of these stars either have such 
low mass loss rates or close companions that their outer 
layers make a disk. Then 0.1 % of all massive star deaths 
would make Type 3 collapsars. For low metallicity stars, 
the latter could be an underestimate. The universal su- 
pernova rate of 6 per second then implies a gamma-ray 
transient event rate of 0.006 per second. If each event 
lasted 10 4 s, there would ~ 50 events in progress at all 
times. Most of these would be at high redshift (low 
metallicity) and the spectrum would be softer and du- 
ration longer than in the lab frame, so the number of 
active sources could be higher. On the other hand, beam- 
ing could increase the brightness and decrease the event 
rate significantly. 

In an extreme limiting case where 10 % of all massive 
star deaths made Type 3 collapsars with duration 10 5 s 
with negligible beaming, one would have ~ 10 5 active 
sources at any one time, or about 1 per square degree. 
If each source radiated 10 51 erg, this could potentially 
contribute to a "diffuse" cosmological background of soft 
gamma-rays. 



The greatest uncertainty in these results is whether the 
matter with high angular momentum will get ejected, by 
mass loss either before the explosion or during the ex- 
plosion, without falling to the center of the star. As 
Table [2] shows, the binding energy of the critical mate- 
rial is very much less than the kinetic energy of a typi- 
cal supernova (~ 10 51 erg), so a weak explosion or even 
strong pulsations during silicon burning could eject it. 
It seems unlikely that this would happen all of the time 
in all sorts of models, so perhaps this enters in chiefly 
as a major uncertainty in the event rates. We also note 
that the converse problem, the optical appearance of a 
"supernova" exploding with only ~ 10 49 erg — 10 50 erg of 
kinetic energy is an interesting one too. One way or an- 
other, it seems unlikely that these stars will all disappear 
without a trace. 
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APPENDIX 
ESTIMATE OF ACCRETION RATE 
Starting from the usual definition of free-fall time, 

r ff = l/^2A^Gp 

we obtain an order-of-magnitude estimate of the accretion rate from 




where we used p(m) = 3m/ (47rr 3 ) and dm = 4Trr 2 pdr, and m is the mass coordinate and r is the radius coordinate 
of a shell. Note that usually p(m) > p(m) in most parts of the star; density inversions can exist at the outer edge of 
convective regions in mixing-length theory especially in RSGs due to recombination and opacity drops. If p <C p, as 
in the case of the outer layers of a BSGs, we obtain M w 2mp/(Tgp) <C M/rg, so BSGs may have a more tenuous 
accretion rate and fainter, softer GRBs than RSGs based on the estimate from rg. On the other hand the average 
density of the star is higher as it is more compact, which may be a slight compensation. Density inversions at the edge 
of RSG envelopes could lead to spikes on accretion rate in this simple approximation. 
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Fig. 1. — Distribution of angular momentum with respect to mass in the presupernova models for V24 and V36. The smooth curves 
show the angular momentum required to form a stable disk at the last stable orbit of a Schwarzschild black hole (lower curve) or Kerr 
black hole (upper curve) including the given mass. The intermediate curve that follows the Schwarzschild curve until far out in the star 
uses the integrated angular momentum in the model to determine the last stable orbit. Where the irregular dark line showing the actual 
angular momentum on the star intersects this line a disk can form. The outer 9 Mq of Model V24 and the outer 10 Mq of Model V36 will 
form a disk. The edge of the helium cores of the two models is apparent in the sharp inflection in the angular momentum at 8.3 Mq and 
15.0 Mq. Mass loss was included in the calculation, but due to the low metallicity, only 0.05 Mq and 0.15 Mq was lost in V24 and V36, 
respectively. Note that, if all the surface material accreted here, the black hole would rotate at nearly its maximum allowed value (i.e., the 
red line intersects the green one at the end). 
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Fig. 2. — Distribution of specific angular momentum in a tidally locked red supergiant in a binary system that has induced corotation 
of the convective envelope. As in Fig. [T] the solid dark line shows the actual distribution of angular momentum while the lighter lines 
show the values required to form a disk around a non-rotating black hole, a Kerr black hole, or a black hole with the integrated angular 
momentum interior to the given mass. The angular velocity in the envelope is constant with u) = 10 — 8 rads - 1 . Note that the black hole 
formed here would be an extreme Kerr hole with a fa 1. 
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Fig. 3. — Final angular momentum for a 250 Mq pair-instability supernova (Model Z250B) in which rotation and magnetic torques were 
included in the calculation. The left panel shows the evolution of the angular momentum and the right panel shows the distribution of 
angular momentum in the presupernova star compared to what is need to make a disk around a black hole. The helium core mass is 166 Mq 
and only matter in the outer 16 Mq has sufficient angular momentum to make a disk. If the outer 16 Mq is not lost during the evolution 
a disk might still form and the black hole could end up rotating rapidly. The radius of this matter is from 1.5 X 10 12 cm to 5.9 X 10 13 cm. 
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Fig. 4. — History of the angular velocity for 8Mq and 16 Mq helium cores in tidally locked binary systems. The two top panels show 
Models 8A (left) and 8B (right) and the bottom panels show Models 16A (left) and 16B (right). Note the tendency of all but the inner 
core to rotate with the (imposed) surface angular speed. At late times the core evolve rapidly and rotates differentially. 



14 





-'6 


I presn 








\ 
N E 




- C depl 








\ 




: C-ign 






\ 


;£ 


17 


- He-depl 












~ Hg burn 








E 










1 


E 










E 


3 






























































1 5 
















/ M He = 8 












cj = 2 x 10* 








"4 















| C depl 




L c-ign 




: 


— — ~* ~- 




w He = 8 : 




cj = 8.5 x 10" J 





18 














- C depl 












I C-ign 








E 


1 7 






















E 












E 


























' 6 












































































M* = 


16 










cj = 2 x 


10" 4 






"4 











5 10 15 

r/M B 



9 - 




Fig. 5. — History of the equatorial angular momentum for 8Mq and 16 Mq helium cores in tidally locked binary systems shown in 
Fig.H 
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Fig. 6. — Final angular momentum for 8Mq and 16 Mq helium cores in tidally locked binary systems compared to that required to form 
a disk around a black hole. Models 8B (upper right) clearly could make a disk and a rapidly rotating black hole. Model 16A (lower left) 
would not. Disk formation in Models 8A and 16B depend on the retention and collapse of a small bit of mass at the surface. 
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Table 1. Summary of Model Characteristics 



Type 


Model 


M disk 


R 


R/v eac 


M 


1 % Mc 2 






(M ) 


(cm) 


(s) 


(10- 4 M Q s- 1 ) 


(lO^crgs- 1 ) 


BSG 


V24 


10.0 


0.4 - 10 X 10 


20,000 


5 


9 




V36 


10.2 


1 - 50 xlO 12 


60,000 


1 


2 




not; 


U. 


l.z X 1U 


1U 






RSG-bin 


S25 


2.8 


0.1 - 9 xlO 13 


3 x 10 6 


0.01 


0.02 


Pair SN 


Z250A 


19.6 


0.9 - 8 xlO 13 


3 x 10 s 


0.7 


1 




Z250B 


20.6 


0.4 - 6 xlO 13 


10 5 


2 


4 




Z250C 


12.1 


0.2 - 2 xlO 13 


3 x 10 4 


4 


8 


WR-bin 


8A 


0.18 


2.5 - 5 xlO 10 


100 


20 


40 




8B 


1.19 


1 - 5 xlO 10 


100 


100 


200 




16A 


0.01 


3 - 4 xlO 10 


100 


1 


2 




16B 


0.83 


2 - 4 xlO 10 


40 


200 


300 



Table 2. Remnant Properties 



Type 


Model 


Period 


-Afpre-SN 


M(j > i cri t) 


BE(j > j crit ) 


a 


A^Hc-corc 








(ms) 


(M ) 


(M Q ) 


(10 49 erg) 




(M ) 




BSG 


V24 


5.4 


23.96 


14.0 - 24 


2.1 


1 


8.32 


0.13 




V36 


3.7 


35.85 


25.7 - 35.9 


1.2 


1 


14.97 


0.11 


RSG-loZ 


025 


5.1 


23.76 






0.52 


10.07 


0.12 


RSG-bin 


R25 


4.6 


12.09 


9.3 - 12.1 


0.0025 


1 


9.19 


0.13 


Pair SN 


250A 




250 


230.4 - 250 


2.9 


1 


142 


0.05 




250B 




250 


229.4 - 250 


8.2 


1 


166 


0.09 




250C 




250 


237.9 - 250 


6.8 


1 


222 


0.13 


WR-bin 


8A 


7.4 


8.043 


7.86 - 8.04 


1.0 




8.04 


0.25 




8B 


3.2 


8.043 


6.85 - 8.04 


7.8 




8.04 


1 




16A 


2.5 


16.087 


16.08 - 16.09 


0.08 




16.08 


0.18 




16B 


1.5 


16.087 


15.26 - 16.09 


12 




16.08 


0.48 



